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Abstract. Fixed income markets share many features with the equity mar- 
kets. However there are significant differences as well and many attempts have 
been done in the past to develop specific tools which describe (and possibly 
forecasts) the behavior of such markets. 

For instance, a correct pricing of fixed income securities with fixed cache 
flows requires the knowledge of the term structure of interest rates. A number 
of techniques have been proposed for estimating and interpreting the term 
structure, yet solid theoretical foundations and a comparative assessment of 
the results produced by these techniques are not available. 

In this paper we define the fundamental concepts with a mathematical 
terminology. Besides that, we report about an extensive set of experiments 
whose scope is to point out the strong and weak points of the most widely 
used approaches in this field. 
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0. Introduction 

Many ideas and techniques developed for the analysis of financial markets can be 
applied both to equities and to fixed income securities. For instance the distribution 
and the dynamical evolution of price variations are studied under the assumption 
that such variations can be described by a random process, regardless of the asset 
category (equity or fixed income) in exam. 

However fixed income markets have many special features, thus specific concepts 
and tools have been introduced in the past to study (and possibly forecast) the 
behavior of such markets. For instance, a correct pricing of fixed income securities 
with fixed cash-flows requires the knowledge of the term structure of interest rates. 

A number of "numerical recipes" have been proposed for estimating and inter- 
preting the term structure of interest rates, yet a solid theoretical foundation of 
such concept and a comparative assessment of the results produced by the different 
techniques is not available. 

In this paper we define formally the concept of absence of arbitrage opportu- 
nities and derive rigorously the existence of discount factors. Moreover the most 
widely used approaches in estimating the term structure are presented and their 
performance is checked with an extensive set of experiments. 

The structure of the paper is the following: section [I] is a reminder of the termi- 
nology used in fixed income markets; section |^ defines a mathematical framework 
for the discount function; section [3] describes the major techniques currently in use 
to estimate the term structure of interest rates; section [| presents the results of our 
numerical experiments; section [| concludes the work with the future perspectives. 



1. Fixed income securities 

A bond is a credit instrument issued by a public institution or a corporate company 
on raising a debt. When the size and the timing of the payments due to the investor 
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are fully specified in advance, the bond falls within the fixed-income securities 
category. 

The main source of uncertainty for a bond markets investor is the default risk 
of the issuer, that is the chance that the issuer becomes unable to fulfill the com- 
mitment of paying all promised cash-flows on time. For this reason the only (true) 
fixed- income securities are those issued by the governments of some developed coun- 
tries or by "well-known" large corporate institutions to which credit rating agencies 
also give their highest rating: these are the borrowers least likely to default on their 
debt repayments. 

Likewise should not be regarded as fixed-income securities those bonds having 
cash-flows indexed to the rate of inflation (index-linked bonds) or those with em- 
bedded optionality giving either the issuer or the holder some discretion to redeem 
early or to convert to another security ^\ . 

The lack of any uncertainty makes fixed-income bonds a useful means for mea- 
suring market interest rates. 

An investor can purchase "new" bonds in the primary market (i.e. directly from 
the issuer) or formerly issued bonds from other investors in the secondary market. 
Furthermore large institutional investors can sell bonds to the issuer as well, that is 
they can pay the price of the bond and undertake to repay the corresponding future 
cash-flows. Hence the bid and ask prices of a bond | are set continuously, as long 
as it remains outstanding, by transactions occurring between market participants. 

In principle a bond can promise any pattern of future cash-flows. However there 
are two main categories. 

Zero-coupon bonds, known also as discount bonds, make a single payment at a 
date in the future known as the maturity date. The size of this payment is called 
the face value, or par value, of the bond. The period of time up to the maturity date 
is the maturity of the bond. US Treasury bills (Treasury obligations with maturity 
at issue of up to twelve months) take this form. 

Coupon bonds make interest payments (often referred to as coupon payments) of 
a given fraction of face value, called the coupon rate, at equally spaced dates up to 
and including the maturity date when the face value is also paid. The frequency at 
which interest payments are made varies from market to market, but generally they 
are made either annually or semi-annually. US Treasury notes and bonds (Treasury 
obligations with maturity at issue of twelve months up to ten years and above ten 
years, respectively) take this form. Coupon payments on US Treasury notes and 
bonds are made every six months. 

In the following, for sake of simplicity, we will focus on US Treasury bond markets 
only. There are a number of advantages in this choice. First of all these markets 
are extremely large regardless of whether size is measured by outstanding or traded 
quantities, are very liquid and have uniform rules. Moreover the number of issues 
is quite high (at present there are about two hundreds obligations outstanding). 
These features make the study of such markets somewhat easier. 

1.1. The Law of One Price. The key idea behind asset valuation in fixed- income 
markets is the absence of arbitrage opportunities, also known as the 'Law of One 
Price'. 



1 To give some examples, we recall the US, UK and Japanese government callable bonds, giving 
the Treasury the option to redeem the bond at face value at any time between two dates specified 
at the time of issue; or the Italian government putable bonds (CTO) giving the holder the option 
to sell the bond back to the Treasury at face value on pre-specified dates; or the UK government 
and corporate convertible bonds giving the holder the option to convert the bond into another 
pre-specified bond at a pre-specified ratio on one or more pre-specified dates (see e.g. jlj]) 

2 respectively the price received when selling and paid when buying that bond 
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The 'Law of One Price' states that two portfolios of fixed-income securities (that 
is two positions taken in the marketplace through buying and selling fixed-income 
securities) that guarantee the investor the same future cash-flows and give him the 
same future liabilities, must sell for the same net price. Phrased another way, a 
portfolio of fixed-income securities giving the investor neither the right to receive 
any future cash-flow nor any future liability, must have zero initial net cost. Con- 
versely a portfolio requiring no initial net investment must give the investor some 
"real" liability, that is some liability whose size is larger than the net amount of 
money earned till then. 

If any violation of these constraints on the prices of all fixed-income securities 
outstanding at a given time occurs, then a profitable and risk-less investment op- 
portunity, called an arbitrage opportunity, or simply an arbitrage, arises. 

Arbitrage opportunities may arise (sporadically) in financial markets but they 
cannot last long. In fact, as soon as an arbitrage becomes known to sufficiently 
many investors, the prices will be affected as they move to take advantage of such 
an opportunity. As a consequence, prices will change and the arbitrage will disap- 
pear. This principle can be stated as follows: in an efficient market [j there are no 
permanent arbitrage opportunities. 

The prices of all fixed-income securities (and in particular those of US govern- 
ment bonds) as they are quoted at any given time, cannot, then, be independent 
one of the other and the existence of (moderate) arbitrage opportunities in financial 
markets can be viewed as a relative mispricing between correlated assets. Such mis- 
pricing can be ascribed to taxation, transaction costs and commissions associated 
with trading that make the net price of a bond different from its quoted gross price 
(and, hence, make the corresponding arbitrage opportunity not really profitable), 
to liquidity effects and other market frictions, or to non-synchronous quotations. 



2. The mathematical setting 

The cash-flows that a US government bond holder is entitled to receive are com- 
pletely determined by the bond face value that, unless otherwise specified, will 
always be supposed equal to $100, by its coupon rate and by its maturity. For this 
reason an outstanding bond can be identified with a triplet (c, m,p), where c > is 
its (semi-annual) coupon rate in percentage points, m > is its maturity in years 
(computed using '30/360' convention) and p > is its gross price in dollars. For a 
US Treasury bill c = and m < 1; for a US Treasury note or bond c > 0. 

We will state, now, formally the condition of absence of arbitrage opportunities 
and derive its consequences on the prices of US government bonds. This will lead 
us to the definition of the discount factors or zero-coupon bond prices. We start by 
giving some definitions. 

Let B := [0, +oo) x R + x R + be the set of all bonds, actually issued by the US 
Treasury or not, and still outstanding at the date in exam. 

Secondly, given a subset S of B, denote by T(S) the vector of the maturities of 
all coupon payments of all bonds belonging to 5, sorted in increasing order. This 
amounts to say that if T(S) = (ti, . . . , t^) for some N € N, then 

{h, . . . , t N } := U( c ,m,p)£s{ m ~ V 2 ; ieNU{0} and m - i/2 > 0} 



3 we recall that according to Malkiel [ |ll[ "A capital market is said to be efficient if it fully and 
correctly reflects all relevant information in determining security prices" 



4 



L. MARANGIO, A. RAMPONI AND M. BERNASCHI 



and ti < ti+i for every i. Moreover for every (c,m,p) 6 S define the cash-flows 
vector ip(c, m) = (ifi(c, m), . . . , (Pn(c, m)) £ M. N by 

!100 + c ,i£m = ti 
c , if ti em - f (i = 1, ... ,iV). 

, otherwise 

Thus (fii(c,m) represents the cash- flow (in dollars) that an investor holding bond 
(c, m,p) is entitled to receive in ti years' time. 

Thirdly, we say that a subset T of B is a complete coupon term structure if and 
only if it contains at least a bond maturing on each coupon payment day (see e.g. 
Pj). Observe that for a complete coupon term structure T the set of the maturities 
of all coupon payments of all bonds belonging to T is the same as the set of the 
maturities of all bonds belonging to T, that is, if T(T) = (tj, . . . , fjv), then 



{ti, ... , t N } :— U( C: „ liP ) Gr {m}. 

Finally we represent a portfolio of B bonds (ci, Tni,pi), . . . , (cb, tubtPb) by a 
vector (qi, . . . , qg) G M B such that for each j, q.j represents the traded quantity of 
bond (cj,rrij,Pj). By definition, qj > means that the investor has bought $100qj 
face value of bond j whereas qj < means that he has sold $100|qj| face value of 
that bond. 

We are, now, in a position to state formally the condition of absence of arbitrage 
opportunities among the bonds belonging to a given set S C B, and to show which 



constraints on their prices such condition implies (sec Theorem 2.5 below) 



Definition 2.1. Let S be a subset of B and let T(S) — (fix,... ,tj^) for some 
N eN. We say thatS satisfies the hypothesis of Absence of Arbitrage Opportunities 
if and only if, given {c\,mi,p\), . . . , [cbtTUbtPb) G S (B £ N) the following three 
conditions are fulfilled: 

(NAi) if (qi, . . . , qs) € R s is such that Y^,f=i Qjficji m j) = 0' then Y^,f=i QjPj = 0/ 
(NA 2 ) if (q%, . . . , qs) 6 R B is such that 



ifi = i 



B 

(J , otherwise 

for some fa > 0, then < J2j=x QjPj < fi! 

B 



(TV A3) if (qi, . . . ,Qb) £ K \{0} is such that Y] ■_ ^ qjpj — 0, then there exists 1 £ 



{1,... ,N} such that J2l=i Sj=i Qjfi( c j -> < °- 

Remark 2.2. The ith component of the vector X^-=i Qjf( c ji m j) represents, if it 
is positive, the total cash- flow that an investor holding portfolio (qi, . . . ,9s) will 
receive in U years' time; if it is negative, the total liability that he will have at that 
date. The quantity ^2f =1 qjPj is the (net) price of the portfolio. 



Hence, condition |(/VAi) means that a portfolio of bonds giving the investor 



neither the right to receive any future cash-flow nor any future liability must have 



zero initial net cost. Conversely (NA3) means that a non-trivial portfolio requiring 
no initial net investment must, sooner or later, give the investor a liability of larger 
size than the total net amount of money received till then. 

The meaning of condition yVAg] is self-evident and will not be discussed any 
further. 
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Definition 2.3. Given B + 1 bonds (c,m,p),(ci, mi, p±), ... ,(cb, tub, Pb) in B, 
a portfolio (31, . . . , qs) of the last B bonds such that X^=i 9iV( c j7 m j) = v( c i m )> 
is called a replicating portfolio of bond (c,m,p). In fact it entitles an investor to 
receive exactly the same cash-flows as if he held bond (c,m,p). If furthermore, 
each bond (cj,rrij,Pj) expires on a different coupon payment date of bond (c,m,p), 
then we will refer to such a portfolio as a minimal replicating portfolio. 



Remark 2.4. Since the case qj = is trivial, condition (NAi) amounts to say that 
the price of each bond (c, m,p) in <S equals the ones of all its replicating portfolios, 
whenever they exist. 

Theorem 2.5. Let T C B be a complete coupon term structu re and let T(T) — 



(t\, . . . ,tff) for some N S N. //Card(T) > N and if condition (NA\) is satisfied, 
then there exist d\, . . . , d/v £ R such that for every (c, m, p) G T 

N 

(2.1) p = y"V^(c, to). 

i=l 

//, furthermore, conditions [N A2~\ and \(N A3) are fulfilled as well, then 1 > d± > 
d 2 > ■ ■ ■ > d N > 0. 

A proof of this theorem is reported in the Appendix. 



Remark 2.6. In order to prove Theorem 2.5 we need to find N bonds in T such that 



the matrix whose columns are their cash-flows vectors have maximum rank N. The 
hypothesis that T forms a complete coupon term structure just helps to ensure that 
this choice is possible. However this is not the only case in which this is true. To 
realize this it is sufficient to think that if S — {(c\, mi,pi), (02, m2,P2), (03, m^,p^)] 
with mi = i, ni2 — TTI3 — I and (of course) C2 7^ C3, then the matrix whose columns 
are the cash- flows vectors of these three bonds has maximum rank 3, even if <S is 
not a complete coupon term structure. 



Remark 2.7. By equation (2.1) with c = it is apparent that each di represents 
the price of a (hypothetical) zero-coupon bond in T with unitary face value and 
maturity m = tj. For this reason d\, . . . , djv are called the discount factors, or 
zero-coupon bond prices, corresponding to the maturities ti,.. . , izv- 

Definition 2.8. Given N discount factors di, . . . , djv corresponding to maturities 
t\, . . . ,tpf, we define the spot rate (of return), or zero-coupon yield, Si, by 

Si := -i Igdi 

ti 

Furthermore, given P S N, if there exists > i such thattj^ — ti — -j? we define 
the P -periods forward rate (of return) Q , ff , by 

h ' P S ' 

Since J- = exp ( f* z Sj dt^ , then Sj represents the continuously compounded rate 
of return per unit time on a ij-years investment, made today, in the bond markets. 
Analogously ff represents the continuously compounded rate of return per unit 
time on a P-periods investment to be made ti years ahead. 

As it is common in financial literature, we will refer to the set of the spot rates 
corresponding to all maturities as the (cross-sectional) term structure of interest 
rates. 



4 as it is common in financial literature, we will refer to any six months' time interval as a 
period 
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2.1. The method of Carleton and Cooper. If a complete coupon term struc- 
ture is available and if conditions (iVAi)| - (NA3) are (strictly) satisfied, then The- 



2.5 tells us how to compute a complete set of decreasing discount factors 



(one for each coupon payment date) that summarize the whole information carried 
about nominal discount rates by bonds in our data set. 



When small violations of condition (NAi ) arc detected in the bond set T forming 
a complete coupon term structure, Carleton and Cooper [^) suggested to consider 
the prices of the bonds in T as random variables with expected values satisfying 
(./VAi)| - [NA3) . According to their hint, equation ( |2.lD must be modified by adding 



a bond-specific error term to (say) its right-hand side: 

N 

(2.2) pj =^2dm(cj,mj) +e h 

i=l 

where j ranges over all bonds in the set T, the e/s are random variables with 
E{ej} ee 0, and 1 > di > d 2 > . . . > dw > 0. 

Estimators, d\, . . . , <i/v, of the discount factors d\, . . . , d/v, subjected to the con- 
dition 1 > di > d- 2 > . . . > dpj > 0, can, then, be attained by a constrained least 
squares procedure with all bonds in the set T . 

We recall, however, that the least-squares estimators d\, . . , , d/v of the d,'s exist 
and are unique provided that the regressors' matrix, <f> = {(<Pi(cj, rrij))), has maxi- 
mum rank N, which is certainly true if T is a complete coupon term structure (see 
Remark ^^). If, on the contrary, the rank of the matrix <I> is less than N, then the 
problem of estimating the discount factors d\ , . . . , d^ by a least squares procedure 
is improperly posed, or ill-posed, in the sense of Hadamard (see, e.g., ||) and does 
not admit any solution. Hence the completeness hypothesis of the term structure 
of T (or a weaker version of its) is necessary to apply the method of Carleton and 
Cooper. 



3. Estimation of the term structure 



In general, the bond set you decide to use for your analysis does not form a 
complete coupon term structure. This is true, for example, for the set we used in 
the present paper. To realize this look at the next figure 1 where we have reported 
the maturity spectrum of all of 206 bonds in our sample. 

Since no bond has maturity of 10 up to 15 years and the maturity of the longest 
(coupon) bond reaches about 30 years, the term structure of coupon bonds in our 
data set cannot be complete. 

In 1971 McCulloch Jl§ introduced a new technique for studying the term struc- 
ture of interest rates based on the following assumption . 



Hypothesis 3.1. If the bond set S C B satisfies conditions (N ~Ai)-(N A3 
there exists a decreasing function d: [0, +00) — > (0, 1] such that d(0) — I, lim t . 
and that for every bond (c, m,p) S S 



then 
+00 d(t) 



N 



P : 



^d(t 4 )^(c,m). 



Remark 3.2. Function d is called the discount function for the bond set S. For 
each t > 0, d(t) represents, then, the present value of $1 to be received in t years' 
time. When S = T forms a complete coupon term structure, the discount factors 

5 Actually McCulloch's statement of Hypothesis differs somewhat from ours: he assumed 
that coupon payments arrive in a continuous stream instead of semiannually. However our formu- 
lation corresponds to the one preferred by most authors after McCulloch 
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Figure 1. 



di, . . . , d/v, whose existence is guaranteed by Theorem 2.5, are just the values that 
d takes on the corresponding maturities t%, . . . ,tpj, that is di = d(ti) for every i. 

The definitions of spot rate and of P-periods forward rate can be easily extended, 
by means of the discount function d, to all maturities t > 0. In particular the spot 
rate function s: (0, +oo) — > K+ is defined by 

(3.1) s(t):=~lgd(t) 



t 



and the P-periods forward rate function f p : [0, +oo) 



by 



f F (t) 



•If 



d(t + P/2) 



P "° d(t) 

Furthermore, if d is differentiable, we can define the instantaneous forward rate 
(of return) function / : [0, +oo) — > R + by 



(3.2) 



d(t). 



Then f(t) represents the simple net rate of return per unit time of an investment 
in the bond markets over an infinitesimal time interval starting t years ahead. 
Moreover the P-periods forward rate f p (t) equals the average of the instantaneous 
forward rate / over the length P/2 interval [t, t + P/2], that is 

r, f t + P/2 

! P {t)^-J t f(s)ds. 
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It is important to understand, however, that the existence of a discount func- 
tion, even when conditions ^NAi) -(NA^) are strictly satisfied, is nothing but a 



hypothesis that scholars and practitioners assume in order to fill the gaps in the 
maturity spectrum of bond samples that do not form a complete coupon term 



structure. It is by no means a consequence of (NAi)-(NA 3 ). As we have already 
pointed out at the end of section 0, building a complete set of discount factors is 
an ill-conditioned problem when the term structure of coupon bonds in the sample 
is incomplete. This explains the difficulties met with by people when trying to 
estimate a discount function from a generic bond data set. 

As to the regularity properties of the discount function, some researchers (see, 
e.g., 0) assert that on the ground of mere economic considerations we can only 
require that d be monotonic decreasing, any further restrictions being not justified. 
On the contrary, most authors assume also that d be twice differentiable so attaining 
a smooth forward rate curve. In particular Langetieg and Smoot JhJ argue that a 
non-smooth forward rate function / should give rise to arbitrage opportunities and 
that, consequently, any irregularity of / should be quickly priced out in an efficient 
market. 



If we assume Hypothesis 3.1, then the problem of estimating the term structure 
of interest rates from a bond data set reduces to choosing an rt-parameter family 
of functions d(-; a.) : [0, +oo) — > R that we consider capable of capturing the main 
characteristics of the bond markets in exam and determining the parameter vector 
a € R ra by minimizing the sum of squared residuals (least squares fitting procedure) 

2 



(3.3) min 



N 



Pj - } y d(tj;a)<pi(cj,mj) 



where j ranges over all bonds (cj,mj,pj) in our data set. The choice of the functions 
d(-; a), that at the very end is always a matter of judgment, is crucial to the quality 
of our fit. 

The test functions d{-; a.) that are usually employed in the fitting procedure ( |3.3| ) 
may be roughly divided in two categories: piecewise (polynomial or exponential) 
functions or spline and functions generated by "parsimonious" models. 

Spline functions are more capable of capturing genuine bends in the discount 
function but have the drawback of over-fitting risk; parsimonious models, on the 
contrary, are less flexible but do a better job in removing noise from the data. So 
the choice of spline functions is better for getting a "good fit" whereas parsimonious 
models should be privileged for sake of "smoothness" . 

3.1. Spline based techniques . The most widely used spline method, specially 
among practitioners, is due to McCulloch []l2| (see also |lj]) and is based on the 
use of polynomial spline functions. 

We recall that, given — oo < K\ < K% < . . . < Kk < +oo (for some k > 2) and 
reNU {0}, we say that / : [K i , Kk] — > R is an r-degree polynomial spline function 
with knot points K\, K2, ■ ■ ■ , Kk, if and only if / is an r-degree polynomial in each 
interval [Ki, Kt+i] and if it is continuous, along with all its derivatives up to order 
r- 1 (if r > 0), in [Ki,K k ]. 

The reasons underlying the choice of polynomial spline functions are as follows. 
It is apparent that a very simple, but somehow naive, approach to the problem of 
fitting a discount function to a bond data sample is to use polynomial test functions. 
However such functions have uniform resolving power, that is they tend to fit better 
at the short end of the maturity spectrum where the greatest concentration of bond 
maturities occurs (see figure Q), and worse at the long end. A possible way out is, 
of course, to increase the order of the polynomial but this can cause instability in 
the parameter estimates. On the other hand, practitioners often require methods 



A CRITICAL REVIEW OF TECHNIQUES FOR TERM STRUCTURE ANALYSIS 9 

endowed with a sufficiently high flexibility to have the chance to choose a priori 
how good their fit of the discount function will be in various regions of the maturity 
spectrum. For instance, they may desire to improve the goodness of their fit in those 
maturity ranges where many bonds are traded. 

A fitting procedure that meets all these requirements involves the use of poly- 
nomial spline functions. By means of a convenient choice of the number and the 
position of the knots, such method allows to achieve the desired resolution degree 
along the whole maturity spectrum, even employing relatively low order polynomi- 
als (which means more stable curves). 

So the choice of the number and the position of knot points affects dramatically 
the quality of the estimated term structure of interest rates. 

McCulloch suggested to choose a number k of knots equal (approximately) 
to the square root of bonds in the sample and to place them in such a way that 
A'i = 0, Kk equals the maturity of the longest bond and that each subinterval 
[Ki, -Ki+i] contains approximately the same number of observed maturities. This 
choice should provide an increasing resolving power as the number of observations 
increases and should avoid over-fitting phenomena. 

Currently, practitioners and financial analysts prefer to place knots at 0, 1, 3, 
5, 7, 11 = 10 + and 30 years, so splitting up the whole maturity range into the 
intervals (0,1), (1,3), (3,5), (5,7), (7,11), and (11,30). The rationale behind this 
choice is that most asset management companies divide their assets into groups (or 
categories) that always lie in one of the above intervals of the maturity spectrum. 
This fact is reflected also in the existence of several benchmark securities with 
maturity of about 1, 3, 5, 7, 10 and 30 years. 

Given k knot points K\, . . . , Kk, the set of all r-degree polynomial spline func- 
tions on [Ki, Kk) form a linear space of dimension k + r — 1. As a matter of fact 
we have the chance of choosing arbitrarily the r + 1 coefficients of the polynomial 
in the first subinterval [K\, K2] and one coefficient for each polynomial in the fc — 2 
remaining subintervals [Ki, Ki + {\ with i = 2, . . . , k — 1. The others parameters are 
fixed by the conditions f^(Ki - 0) = f {n) {K t + 0) for all i e {2, . . . , fe — 1} and 
aline {0,... ,r-l}. 

Let {/o,/i, ■ • • , /r+fc-2} be a base for the linear space of r-degree polynomial 
spline functions with knot points K\, . . . , Kk- Without loss of generality we can 
assume that fo(t) = 1 (as a constant function is a polynomial spline function) and 
that /i(0) = . . . = / r +fc-2(0) = 0. Then a = (a%, . . . , av+fc-2) and the McCulloch 
model for the discount function is 

r+fc-2 

d(t;a) = l+ Y, "j/jW- 

i=i 

The corresponding least squares problem is therefore linear in the parameter vector 
a. 

A major limitation of this method is that it does not force a monotonous de- 
creasing behavior of the resulting discount function. As a consequence, meaningless 
results, like negative forward rates, may be obtained. In this respect Barzanti and 
Corradi || (see also || and Q]) have recently, proposed the introduction a set of 
linear constraints on the parameters a±, . . . , a r +k-2 that ensure the desired mono- 
tonicity of the employed spline functions. 

A further criticism directed by Vasicek and Fong |l^] to McCulloch's method is 
that most general-equilibrium models of the term structure of interest rates (see, 
e.g., ]l8| and 0) expect an exponential form for the discount function. However 
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piecewise polynomials have a different curvature compared to an exponential func- 
tion. Hence they conclude that "conventional" polynomial splines cannot provide 
a good fit to an exponential-like discount function. 

Following the arguments of Vasicek and Fong, Langetieg and Smoot jlO| and 
Coleman, Fisher and Ibbotson Q proposed, respectively, to fit a polynomial spline 
to the spot rate function and to the instantaneous forward rate function rather 
than to the discount function ^] . They argued that, assuming an exponential form 
for the discount function, their regression procedure should give better results (i.e. 
more "realistic" from the financial viewpoint) compared to McCulloch's technique. 

By ([H]) and |T2|) d(t) = exp(-is(t)) and d(t) = exp (- f* f(s) ds^j . Hence the 



method by Langetieg and Smoot amounts to take a = (ao, . . . , av+fc-2) (for some 
reNU {0} and some k > 2) and 

(r+k-2 

where {/o,/i,... , f r+k-2} is a base for the linear space of r-degree polynomial 
spline functions with knot points K\ 1 . . . , Kk- 

Analogously, since the integral of an r-degrce polynomial spline function is an 
(r + l)-degree polynomial spline, the method by Coleman, Fisher and Ibbotson 
amounts to take a. — (ao, . . . , a r +k-i) and 

(r+fc— 1 

where, as above, {/o, /1, ■ ■ ■ , /r+fc-i} is a base for the linear space of (r + l)-dcgrcc 
polynomial spline functions with knot points K\ , . . . , Kk and with the property 
that f (t) = 1 and that / x (0) = . . . = / r +*-i(0) = 0. 

The methods by Langetieg and Smoot and by Coleman, Fisher and Ibbotson 
share the same limitation of the method by McCulloch in that it is impossible 
to force a monotonous decreasing behavior of the resulting discount function by 
simply introducing linear constraints on the parameters to be determined. The 
only exception is for the method by Coleman, Fisher and Ibbotson with r = 0, that 
is when a piecewise constant functions are used to fit the instantaneous forward 
rate function. Hence, again, financially unacceptable results may be found. 

3.2. Parsimonious models. It is readily apparent that, although spline based 
techniques are still widely used to estimate the term structure of interest rates, 



they present many drawbacks (as pointed out, for example, by Shea 15, |T6|]). As 
a matter of fact several choices, such as the number and the position of knots and 
the degree of the basic polynomials, affect dramatically the quality of fit. 

An alternative approach to the estimation of the term structure of interest rates 
is the use of the so-called parsimonious models. The main idea behind this class of 
models is to postulate a unique functional form on the whole range of maturities 
for the discount function. This implies that in parsimonious models some basic 
properties of the discount function, such as monotonicity, we are interested in can 
be imposed a priori. Furthermore, the number of parameters to estimate is typically 
much smaller than that needed for a typical spline based approach. 

Let us describe briefly some of the parsimonious models proposed in the litera- 
ture. 



"Actually Coleman, Fisher and Ibbotson proposed only the use of piecewise constant functions, 
that is of 0-degree polynomial splines. However their method can be easily extended to polynomial 
splines of any order. 
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In M the functional form for the discount function is proposed to be 



d(t;a) = exp(— ^""^ a.jt 3 ). 

3=1 

The unknown parameters a\, . . . ,a n can be estimated either by (non-linear) 
least square regression or by a maximum likelihood approach. The latter one is 
found to perform better, probably, because of the explicit consideration of the 
heteroskedasticity in the data. Notice that this model is a particular case of the 
ones by Langetieg and Smoot and by Coleman, Fisher and Ibbotson with just two 
knots. 

A particularly relevant approach is the one proposed by Nelson and Siegel [ fl4| . 
They attempted to model explicitly the implied instantaneous forward rate function 
/. The Expectation Hypothesis provides a heuristic motivation for their model since 
/ is modeled as the solution of a second order differential equation. The functional 
form suggested for / is 

(3.4) f(t; a) = p + Pie^ + P^-e^, 

T 

where a = (/3q, Pi, /?2 , t). It implies the following model for the discount function: 
d(t;a) =exp (-t (#, + (ft + ) y (l - e^ T ) - p 2 e~ t/T ) ) . 

The flexibility of this model is explained by the authors by observing that the 
instantaneous forward rate is the linear combination of three components, f s (t) :— 
e - */ 7 ', f m (t) '■= Te~ t l T and fi(t) := 1, modeling respectively the short-, medium-, 
and long-term behavior of /. The parameters Po, Pi and P2 measure respectively 
the strength of each component, whereas r is a time-scale factor. 

In jlTj an extended model for the instantaneous forward rate / was proposed to 
increase the flexibility of the original model by Nelson and Siegel model. It was 



obtained by adding a fourth term in equation (3.4), with two extra parameters, 



which allow the instantaneous forward rate curve for a second "hump" : 
f(f, a )=po+P 1 e- t ^+p 2 -e- t ^+P 3 -e- t / T \ 

n t 2 

The model proposed by a J. P. Morgan group and reported by Wiseman [ po| , 
is the so-called Exponential Model. The instantaneous forward rate function is 
modeled by 

n 

f(t; a) a j exp(-6^), 

3=0 

where a%,bi, . . . ,a n ,b n are the parameters to be estimated. It has been noticed 
that this model is able to capture the macro shape of the forward rate curve rather 
than the features of a single data set. 

4. Numerical Experiments 

In this section we report the results of a set of numerical experiments for estimating 
the term structure of interest rates from a cross-sectional data set of US government 
bonds. 

The data set was made available by Datastream. It contains the annual coupon 
rate, the maturity date and the gross (or dirty or tel quel) price for all (fixed- 
income) US Treasury bills, notes and bonds outstanding at the dates from May 
31 st 1999 to June 11 th 1999 (two weeks). 

7 i.e. the sum of the quoted, or clean, price and the accrued interest (computed by the '30/360' 
convention) 
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Figure 2. 



Hereafter, for economy of presentation, since the results are very similar, only the 
results for June 3 rd 1999 (when 33 bills and 173 notes and bonds were outstanding) 
will be presented and analyzed. 

As first step, we built all minimal rep licating portfolios (if any) of each available 
bond to double check whether condition (NAi) is actually violated by bonds in the 
data set. As a "measure" of the arbitrage opportunities for a bond, we took the 
largest absolute value of the difference between the price of the bond and the price 
of all its minimal replicating portfolios. The results of this test are reported in the 
figure ||. 

Note that just 69 bonds, out of 206, admit a (minimal) replicating portfolio, the 
maturity of such bonds reaching at most 7 years. This is easily seen by comparing 
figures H and [|. Moreover arbitrage opportunities amount to less than $0.1 for 
bonds with maturity up to 3 years, and range between $0.2 and $1 for bonds with 
maturity within the range 3 — 7 years. 

To estimate the term structure of interest rates we have resorted to the ap- 
proaches suggested by Carleton and Cooper ||, to the spline approaches by Mc- 
Culloch [jl2| [[3] and by Coleman, Fisher and Ibbotson |7| and to the parsimonious 
model by Nelson and Siegel [[ti] (see section || for a description of these methods) . 
The technique proposed by Svensson |Q gave us the same kind of results as that 
proposed by Nelson and Siegel. Hence this set of experiments will not be reported 
in the present paper. 
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To apply the method of Carleton and Cooper to a bond sample, <S, whose ma- 
turity spectrum contains "gaps", it is necessary to shrink the sample itself to a 
suitable subset, T, that forms a complete coupon term structure. After that we 
can estimate the discount factors by a constrained least squares procedure as de- 
scribed in section ||. In our case the subset T contains 152 (out of 173) US Treasury 
notes and bonds (the longest of which has maturity of about 7 years) and, obviously, 
all of 33 US Treasury bills. 

For our numerical experiments, neither the covariance matrix, cry = E{si£j} nor 
other information about the statistical properties of the errors ej's was available. 
For such reason we have resorted to a straight least squares procedure, instead 
of a weighted one, so as to obtain unbiased, even though not minimum variance, 
estimators of the discount factors. Our results are reported in figure^. 

Note that the discount factors approach one as their maturity tends to zero. The 
discount factors vary quite regularly with maturity, unlike spot rates and, even more 
significantly, one-period forward rates. These present irregularities for ma turities 



of about 1 year and above 3 years, where the largest violations of condition (NAi) 
occur. The plot of the residuals 

N 

Pj -^rdi<Pi(cj,mj), 
i=l 

shows that, as expected, the fitting procedure performs worse where condition 
(jVAjj is more severely violated. 



McCulloch's method has been tested on our data sample with two different sets 
of knot points. The first set (hereafter referred to as "knots ala McCulloch" ) is 
built following McCulloch's indication and consists of 14 knots placed at 0, 0.2, 0.4, 
0.6, 0.9, 1.4, 1.8, 2.4, 3.4, 4.3, 6.3, 16.6, 22.8, 29.7 years. The second set is made of 
seven knots placed at 0, 1, 3, 5, 7, 11 and 30 years. The base we used (see section 



3.1) consists of the following functions: 
/«(*)=*' (* = 0,...,r-l) 1 

[0 , ift<Jf i _ r+1 , 

fi(t) = <(t- K^ r+l ) r , if Ki-r+t < t < Ki-r+2, 

[(t - Ki- r+1 y - (t - Ki- r+2 y , ifK^ r+2 <t 

(i = r, . . . , r + k — 2). 

and polynomial spline functions up to the fifth degree have been used. Our results 
are reported in figure || for the first set of knots, and in figure || for the second set. 

These graphs show all drawbacks of McCulloch's method. First the outcome of 
the regression procedure depends strongly on the positions of the knots and on the 
degree of the polynomial spline employed. This feature is apparent in the forward 
rate curves (compare figure ^ and whereas the fits for the discount function and 
the spot rate are more stable. Moreover, a polynomial spline is not, in general, 
a decreasing function, hence negative forward rates may occur. It is interesting 
to note how results for one-period forward rates obtained both by the method 
of McCulloch (or Coleman, Fisher and Ibbotson as we will see soon) and by the 
method of Carleton and Cooper share the same oscillating behaviour within 3 and 
7 years. This suggests that such behaviour actually depends on the violations of 
the condition of absence of arbitrage opportunities that occur at those maturities 
(see figure ||) . 

To assess the criticism of Vasicek and Fong [|l9] to McCulloch's technique (see 
section [|) , we assumed an exponential discount function of the form 

(4.1) d(t) = e- a06t , 
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which corresponds to a spot rate function and a forward rate function identically 
equal to 6% per year, and generated two sets of artificial data as follows. For each 
bond (cj,mj,pj) in the original data set, S, we defined a fake price: 

N 

where T(S) = (ti,... , ijv) and ej = for the first set of artificial data (named 
hereafter "exact" data) whereas Ej ~ A/"(0, 1) for the second set ("noisy" data). 

After that McCulloch's method, with knots ala McCulloch, has been applied to 
both "exact" and "noi sy" bond data so attaining an estimate, d, of the discount 
function d defined in (4.1). The results reported in figures ^ and (?] seem to con- 
firm Vasicek and Fong's opinion. McCulloch's method performs pretty well on the 
"exact" data, whereas it provides a bad fit when "noisy" data are used, even if the 
artificial perturbations introduced are quite small compared to the typical price of 
a bond (< 1% on average). 

The method of Coleman, Fisher and Ibbotson has been tested with the same 
sets of knot points and the same base of spline functions that we used to test 
McCulloch's method. Polynomial spline functions of various degree have been used. 
For 0-degree spline functions, that is when a piecewise constant function is fitted 
to the instantaneous forward rate function, the least squares procedure has been 
constrained so as to attain a monotonic decreasing discount function. Our results 
are reported in figure || for knots ala McCulloch, and in figure [)] for knots placed 
at 0, 1, 3, 5, 7, 11 and 30 years. 

These graphs show that the method of Coleman, Fisher and Ibbotson shares 
the same drawbacks of McCulloch's method. However the risk of getting negative 
forward rates seem to be reduced. 

The graphs in figure [To] are produced following the approach of Nelson and Siegel. 
The main features of the estimates based on parsimonious models are readily seen. 
Although for the discount function the differences with other approaches are pretty 
mild, the constraints imposed on the functional form of the instantaneous forward 
rate function get rid of the oscillating behavior peculiar to spline based techniques. 

We conclude this section by showing a table where the Root Mean Square Error 
(RMSE) obtained in our numerical experiments is reported. This allows us to 
compare the quality of the fit among all employed methods (for various choices of 
the parameters, such as knot points or degree). 

As expected the quality of the fit is higher for spline based techniques and lower 
for parsimonious models. Moreover it is worth noting that it is not true that 
increasing the order of the spline or the number of knots yields a better fit. 



5. Conclusions 

We have presented a critical review of the main methods for the estimation of the 
term structure of interest rates in fixed income markets. 

It is apparent that most of the activity in this field lacks of a solid mathematical 
basis and there is not a single technique which is definitely better than others. 

A number of elements contribute to form this scenario: 

Bad practice: there is a general trend of re-cycling existing techniques without any 
critical evaluation of their real effectiveness. 

Low-quality data: errors can be present in the input data (wrong prices). A large er- 
ror is immediately detected whereas smaller errors may induce to take into account 
"fake" arbitrage opportunities. 



8 i.e. the square root of the arithmetic mean of squared residuals, see equation (3.3) 
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Figure 3. 
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Figure 4. 
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Method 


Knot points 


spline degree 


RMSE 


Carleton 






.0807 


McCulloch 


ala McCulloch 


2 


.2428 


McCulloch 


ala McCulloch 


3 


.2441 


McCulloch 


ala McCulloch 


4 


.1896 


McCulloch 


0, 1, 3, 5, 7, 11, 30 


2 


.2493 


McCulloch 


0, 1, 3, 5, 7, 11, 30 


3 


.2358 


McCulloch 


0, 1, 3, 5, 7, 11, 30 


4 


.2389 


McCulloch 


0, 1, 3, 5, 7, 11, 30 


5 


.1846 


Coleman 


ala McCulloch 





.3003 


Coleman 


ala McCulloch 


1 


.2350 


Coleman 


ala McCulloch 


2 


.2328 


Coleman 


ala McCulloch 


3 


.1926 


Coleman 


0, 1, 3, 5, 7, 11, 30 





.4560 


Coleman 


0, 1, 3, 5, 7, 11, 30 


1 


.2242 


Coleman 


0, 1, 3, 5, 7, 11, 30 


2 


.2405 


Coleman 


0, 1, 3, 5, 7, 11, 30 


3 


.2282 


Nelson 






.5102 



Table 1. Fit quality 



Differences among the markets: fixed income markets are not homogeneous. The 
number of data available for a cross-section analysis of the U.S. Treasury securities 
is very different, for instance, from the number of outstanding Italian government 
bonds. The features of a market should be considered very carefully before applying 
any technique. 

We expect to extend our work in two directions: 
The development of algorithms to find the best position of the knots for spline- 
based interpolation methods. We consider genetic algorithms very promising for 
this kind of problem. 

A mathematical analysis of equilibrium term structure models. This is pretty im- 
portant since there is already evidence of relevant differences between the predic- 
tions of existing models and the experimental data (an example is the volatility 
term structure). 

The importance of "intuition" in financial markets analysis cannot be under- 
estimated. However a more precise distinction among hypotheses, mathematical 
implications of these hypotheses and empirical observations is definitely required 
in this field. 

6. Appendix 



This appendix is entirely devoted to the proof of Theorem 2.5 that we write 
again here for our convenience. 

Theorem 6.1. Let T C B be a complete coupon term structure and let T(T) — 
(ti, . . . ,t]\[) for some iVeN. //Card(T) > N and if the following condition: 
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Figure 6. 



ifi = i 
otherwise 



(NAi) if(qi, . . . , gs) £ R B are such that J2f=i Qj l f(. c ji m j) = 0? then J2f=i IjPj = ®> 
is satisfied, then there exist d\, . . . , aV £l such that for every (c, m,p) € T 

iv 

(6.1) p =y]diipi(c,m). 

i=l 

If, furthermore, conditions: 
{NAi) if (<?i, . . . , gs) 6 K B is smc/i i/iai 

3=1 I 

/or some > 0, then < X)f=i < ft! 
{NA3) if {qi,... ,qs) £ R s \{0} are swc/i that ^2f = iqjPj = 0, then there exists 
% £ {1, . . . ,iV} smc/i that J2l=i I2j=i Qj<Pi{ c ji m j) > °; 
are fulfilled as well, then 1 > di > 0*2 > . . . > aV > 0. 

Proof. Let (ci, mi,pi), . . . , (cjv, tun,pn) £ T be such that rrij = for every j g 
{1, . . . , N}. The existence of A*" bonds like that in the set T is ensured by the 
hypotheses that Card(T) > N and that T forms a complete coupon term structure. 
Let $ = be an N x TV matrix such that 

(6.2) $y := <pi(cj,mj) 
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Figure 7. 



for alH,j G {1,... , N}. 

Finally let (c, m,p) be an arbitrary bond of T. 

Since $ is an upper triangular matrix with no vanishing clement on its principal 
diagonal, it is non-singular. 

Now let qi, . . . , q^ S K be such that 

N 

3»=S( $_:l )iiVi(c,m) (i = l,...,iV) 



i=l 



Then 



^^^(cj.mj) = <p(c,m). 



By condition (ATAi), this implies equation (6.1) with 



iV 

di:=X)Pi( $ "% (i = l,...,JV). 
i=i 

In order to conclude the proof of the first part of the theorem, we have to show that 
the di's are independent of the choice of (ci,mi,pi), . . . , (cat, itin,pn) such that 
rrij = tj for every j if several such choices are possible. To do this, let j £ {1, . . . , N} 
and let (cj,rhj,p~j) £ 7~\{(cj, mj, Pj)} be such that TOj = ij. 
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Discount function (by method ot Langetieg & Smoot) 
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Figure 10. 
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(6.3) 



Let $ = (($ij)) be an N x N matrix such that 

<Pi{cj,mj) , if J 7^ J 

for all z, j € {1, . . . , N}, and let 

JV 

di := + (i = 1,... ,JV). 

5=1 

i/j 

We want to show that di = di for every i G {1, . . . , iV}. Since both J and rhj,pj) 
has been arbitrarily chosen, this will conclude the proof of the first part of the 
theorem. 



Observe that $.j = for all j ' ^ j and that by equations (6.1), ( |6.2| ) and ( |6.3| ) 



A' 



Pi = 



(i = i,...,AT) 



and 



A 



Pi 



E d ^r 



Then for every i G {1, . . . , iV} 

AT 

3=1 

A 

= E* 



A 



hi 



AT 



E^i^" 1 ) 

i=i 



= di 



In order to prove the second part of the theorem, let (ci, mi,pi), . . . , (cjv, itin,pn) € 
T be such that mj = i,- for every j £ {1, . . . , N} and let (gi, . . . , gjv) E be a 
portfolio of these bonds such that 



A 



ft , if i = t 
E*^( c i' m i) = { &+ 1 > if * = « + 1 



i=i 







otherwise 



for some i e {1, 
is 



, N} and some /j+i G R. By (6.1), the price of this portfolio 



A r 



i=i 



If /j > and = 0, then by condition {NA2), < d^fi < /j, which implies that 
< di < 1 and, since z has been chosen arbitrarily, that 

< d, < 1 (i = l,... ,N). 

Now suppose, by contradiction, that di < dj+i and let fi, /j+i be such that 

N 

difi + dr+iA+i = gjPj = 0. 
i=i 

Since a\, dj+i > 0, we can assume, without loss of generality, that fi > and that 
fi + fi+i > 0. However condition (TV A3) implies that either /, < or /» + /»+! < 0. 
This is the wanted contradiction. □ 



a critical review of techniques for term structure analysis 25 

References 

[1] Anderson, N., F. Breedon, M. Deacon, A. Derry, and G. Murphy, Estimating and Interpreting 

the Yield Curve, John Wiley and Sons, Series in Financial Economics and Quantitative 

Analysis, Chichester, New York, Brisbane, Toronto and Singapore, 1997. 
[2] Barzanti, L., and C. Corradi, A note on direct term structure estimation using monotonic 

splines, submitted for publication to A.M.A.S.E.S. 
[3] , A note on interest rate term structure estimation using tension splines, Insurance: 

Mathematics and Economics 22 (1998), 139-143. 
[4] Campbell, J. Y., A. W. Lo, and A. C. MacKinlay, The econometrics of financial markets, 

Princeton University Press, Princeton, New Jersey, 1997. 
[5] Carleton, W. T., and I. A. Cooper, Estimation and uses of the term structure of interest 

rates, The Journal of Finance 31 (1976), no. 4, 1067-1083. 
[6] Chambers, D. R., W. T. Carleton, and D. W. Waldman, A new approach to estimation of the 

term structure of interest rates, Journal of Financial and Quantitative Analysis 19 (1984), 

no. 3, 233-252. 

[7] Coleman, T.S., L. Fisher, and R. G. Ibbotson, Estimating the term structure of interest rates 
from data that include the prices of coupon bonds, The Journal of Fixed Income (1992), 
85-116. 

[8] Cox, J. C, J. E. Ingersoll, and S. A. Ross, A theory of the term structure of interest rates, 

Economctrica 53 (1985), 385-408. 
[9] A. Kirsch, An introduction to the mathematical theory of inverse problems, Applied Mathe- 
matical Sciences, vol. 120, Springer- Verlag, New York, 1996. 

[10] T. C. Langetieg, and J. S. Smoot, Estimation of the term structure of interest rates, Research 
in Financial Services 1 (1989), 181-222. 

[11] B. Malkiel, Efficient Market Hypothesis, in Newman, P., M. Milgate, and J. Eatwell (eds.) 
New Palgrave Dictionary of Money and Finance, Macmillan, London. 

[12] J. H. McCulloch, Measuring the term structure of interest rates, The Journal of Business 44 
(1971), 19-31. 

[13] , The tax-adjusted yield curve, The Journal of Finance 30 (1975), no. 3, 811—830. 

[14] C. R. Nelson, and A. F. Siegel, Parsimonious modeling of yield curves, Journal of Business 

60 (1987), no. 4, 473-489. 
[15] G. S. Shea, Pitfalls in smoothing interest rate term structure data: Equilibrium models and 

spline approximations, Journal of Financial and Quantitative Analysis 19 (1984), no. 3, 253— 

269. 

[16] , Interest rate term structure estimation with exponential splines: A note, The Jour- 

nalof Finance 40 (1985), no. 1, 319-325. 

[17] L. Svensson, Estimating and interpreting forward interest rates: Sweden 1992-94, Interna- 
tional Monetary Fund Working Paper (1994), no. 114. 

[18] O. A. Vasicek, An equilibrium characterization of the term structure, Journal of Financial 
Economics 5 (1977), 177-188. 

[19] O. A. Vasicek, and H. G. Fong, Term structure modeling using exponential splines, The 
Journal of Finance 37 (1982), no. 2, 339-348. 

[20] J. Wiseman, European fixed-income research: Second edition technical specification - the 
exponential yield curve model, JP Morgan (1994). 

ISTITUTO PER LE APPLICAZIONI DEL CALCOLO "MAURO PlCONE" - C.N.R. - VlA DEL POLICLIN- 

ico, 137 - 00161 Rome - ITALY 



